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ABSTRACT 

We present the first detection and mapping of the HD 32297 debris disk at 1.3 mm with the 
Combined Array for Research in Millimeter- wave Astronomy (CARMA). With a sub-arcsecond beam, 
this detection represents the highest angular resolution (sub)mm debris disk observation made to 
date. Our model fits to the spectral energy distribution from the CARMA flux and new Spitzer MIPS 
photometry support the earlier suggestion that at least two, possibly three, distinct grain populations 
are traced by the current data. The observed millimeter map shows an asymmetry between the 
northeast and southwest disk lobes, suggesting large grains may be trapped in resonance with an 
unseen exoplanet. Alternatively, the observed morphology could result from the recent breakup of a 
massive planetesimal. A similar-scale asymmetry is also observed in scattered light but not in the 
mid-infrared. This contrast between asymmetry at short and long wavelengths and symmetry at 
intermediate wavelengths is in qualitative agreement with predictions of resonant debris disk models. 
With resolved observations in several bands spanning over three decades in wavelength, HD 32297 
provides a unique testbed for theories of grain and planetary dynamics, and could potentially provide 
strong multi-wavelength evidence for an exoplanetary system. 

Subject headings: circumstellar matter — planetary systems: formation — planetary systems: proto- 
planetary disks — stars: individual (HD 32297) 



1. INTRODUCTION 

Debris disks provide the principal means of studying 
the formation and evolution of planetary systems on 
timescales of 10—100 Myr. Evidence for exoplanets in 
these systems can be found by matching density varia- 
tions in debris disks to theoretical models of the gravi- 
tational perturbations caused by planets (e.g., Reche et 
al. 2008). A modest sample of debris disks have now 
been imaged in the visible, and many show substruc- 
ture such as clumps, warps, and offsets, consistent with 
dynamical perturbations by massive planets. However, 
a wide variety of other mechanisms can produce simi- 
lar structures (Moro-Martin et al. 2007, and references 
therein). As different wavebands are sensitive to differ- 
ent grain sizes, which are in turn subject to different 
dynamical influences, multi-wavelength observations of- 
fer the most promising path towards definitively classi- 
fying the physical mechanisms at work in these systems. 
At present, only a few debris disks have resolved ob- 
servations spanning more than a decade in wavelength. 
A particularly critical, though technologically challeng- 
ing, deficit of observations lies at (sub)millimeter wave- 
lengths, which trace large grains primarily affected by 
gravitational forces. To date, bolometer arrays have re- 
solved 8 debris disks in the (sub)mm (e.g., Holland et 
al. 1998, Greaves et al. 1998). However, such low- 
resolution (#beam ^ 10") single-dish measurements are 
limited to the largest, nearest disks. Higher resolution in- 
terferometric observations are needed to access the larger 
debris disk population already imaged at shorter wave- 
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lengths. Some pioneering work has been done in this 
area; OVRO and PdBI have detected and resolved two 
debris disks (Vega: Koerner et al. 2001, Wilner et al. 
2002; HD 107146: Carpenter et al. 2005). Recently, 
ICorder et all (|2008h mapped HD 107146 with the Com- 
bined Array for Research in Millimcter-wave Astronomy 
(CARMA) at 1.3 mm, providing the highest fidelity in- 
terferometric debris disk map to date. Here, we report 
the near-simultaneous CARMA detection of HD 32297, 
the third debris disk mapped with a (sub)mm interfer- 
ometer. With a sub-arcsecond beam, this detection is 
the highest angular resolution (sub)mm debris disk ob- 
servation made to date. 

HD 32297 is a ~ 30 Myr A-star at 112±^ pc 
(|Perrvman et a l. 1997), first discovered to host a re- 
solved debris d isk with HST/NICM OS near-infrared 
(NIR) imaging (jSchneider et al.ll2005h . The discovery 
image showed an edge-on debris disk extending to 400 
AU (3.3"), with an inne r-disk b rightn ess asymmetry in- 
ward of 60 AU (0.5"). iKalasI (|2005l ) subsequently im- 
aged HD 32297 in the optical, revealing an asymmet- 
ric, extended outer disk (^1700 AU, 15" ) likely in t eract - 
ing with the interstellar medium. Later, Rcdficld (|2007l ) 
detected circumstellar gas in this system, reporting the 
strongest Na I absorption measured toward any known 
debris disk. M ost recently, [F itzgerald et al. (2007), here- 
after F07, and lMoerchen et all (|2007| ) resolved HD 32297 
in mid-infrared (MIR), thermal emission. Detailed anal- 
ysis of the spectral energy distribution (SED) by F07 
showed that multiple grain populations may be present 
in the disk. The lack of long wavelength data needed 
to characterize the large grain properties of HD 32297 
motivated the CARMA observations presented here. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed HD 32297 with CARMA, located 7200 
ft above sea level outside Big Pine, California and con- 
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sisting of six 10.4-m and nine 6.1-m antennas, previ- 
ously comprising the OVRO and BIMA arrays. We 
used CARMA's D configuration (baselines: 11—148 m) 
on 2008 March 06 and the more extended C configu- 
ration (baselines: 26—370 m) on 2007 October 26 and 
2007 November 08. We tuned the receivers to a cen- 
tral frequency of 227 GHz. The total continuum band- 
width is 1.5 GHz, contained in three 500 MHz bands. 
Our total observation time was 19 hours in good weather 
with rms path errors of < 175/Ltm and zenith opacities of 
T230 < 0.35. Throughout our observations, we used opti- 
cal offset guiding to maximize pointing accuracy (Corder, 
Carpenter & Wright 2008, in prep). 

We calibrated the data using the MIRIAD software 
package. We performed passband calibration with 15 
minutes integrations on a bright quasar (3C 84, 3C 111) 
observed at the start of each track. We derived time- 
dependent phase solutions from 3 minutes integrations on 
J0530+135, observed following each 15 minutes integra- 
tion on source. In addition, we used 1 minutes integra- 
tions on a secondary calibrator, 3C 120 (observed every 
cycle), to test the astrometric accuracy of the interfer- 
ometer and the integrity of our phase solutions. Finally, 
we flux-calibrated the data using a 5 minutes integration 
on a planet (Uranus, Mars) observed once per track. The 
systematic uncertainty in CARMA's absolute flux scale 
is ^20% (W. Kwon, private communication). 

The 1.3 mm map of HD 32297 combining all obser- 
vations is shown in Figure [T] For this map, we adopted 
natural weighting to provide optimal sensitivity and pro- 
cessed the visibility data as a mosaic to accomodate the 
heterogeneous array. We deconvolved the di rty map us- 
ing th e Steer CLEAN algorithm for mosaics ([Steer et al.l 
1984), set to bring the rms in the cleaned region of the 
residual map to that measured in an off-source region. 
To emphasize possible resolved structure, we restored 
the map with a circular beam with radius equal to that 
of the semimajor axis of the naturally weighted beam 

(Abeam = 0.9"). 

In addition to the new CARMA data presented here, 
we also extracted Spitzcr MIPS photometry from direc- 
tor's discretionary time program 225 (PI: G. Schneider). 
Using the MOPEX package, we manually calibrated the 
24 /an data to remove the strong background gradient 
and "jailbars" evident in the pipeline processed image. 
The 70 jum data did not suffer severe artifacts; in this 
case, we simply mosaicked the data. The 160 /im data are 
contaminated by a spectral leak, which occurs for bright 
sources and is typically corrected using observations of a 
calibration star. However, no such calibration observa- 
tion was made in this program; we, therefore, include the 
160 jtzm data in this analysis only as an upper-limit. The 
MIPS fluxes as derived from aperture photometry are 
#,(23.68 /im) = 0.21 ± 0.01 Jy, #,(71.42 fim) = 0.85 
± 0.06 Jy, #(155.9 /mi) < 0.46 ± 0.06 Jy. In the SED 
modeling (§4), we adopt calibration uncertainties listed 
in the MIPS data handbook and apply a color correction 
appropriate to the modeled grain temepratures. 

3. RESULTS 

Figure [1] shows a 1.3 mm continuum source detected at 
the 7cr level. The observed peak in the map is 2.8 ± 0.4 
mjy beam -1 . The probability of detecting an unrelated 
background source within 1" of the HD 32297 stellar po- 
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Fig. 1.— CARMA 1.3-mm continuum map of HD 32297. The 
black cross marks the stellar position, with the full-width repre- 
senting four-times the total rms millimeter positional uncertainty. 
Contours begin at 2<r and increase by Icr thereafter (<r = 0.44 mjy 
beam -1 , 6>beam = 0.9"). The source morphology suggests that HD 
32297 may be marginally resolved. The measured position angle 
is consistent with that observed in the NIR and MIR (45 — 57°). 
The centroid of the disk emission is offset from the stellar position 
at the 4cr level. 
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Fig. 2. — Vector-averaged visibility amplitudes as a function of 
projected baseline length from phase center. The combined data 
amplitudes decrease with distance from phase center, suggesting 
the source is resolved at the 2a level. The D-array data, on the 
other hand, are consistent with a point-source model (grey dashed 
line), allowing a robust flux measurement. 

sition is <Cl%, based on recent source counts at 850 /im 
(e.g., Scott et al. 2002). Thus, we conclude that the de- 
tected emission is associated with the HD 32297 debris 
disk. The map morphology suggests the source may be 
resolved (see below). However, in the D-array data set 
alone the visibility amplitudes are constant with baseline 
length (Figure [2]), suggesting the source is unresolved. 
We therefore, adopt an integrated flux measured using 
the D-array data set only: 5.1 ± 1.1 mjy. 

The source centroid in Figure [1] appears offset from the 
Hipparcos stellar position (located at phase center). To 
quantify this result, we fit a point source model to the 
visibility data; the measured offset is Ar = 0.43 ± 0.08" 
(Ac* = -0.30±0.08", AS = -0.31±0.07"). This offset is 
significantly larger than that observed for test calibrator, 
3C 120, observed each source cycle (Ar = 0.041 ±0.001", 
Aa = -0.010 ± 0.002", A<5 = 0.040 ± 0.001"). In ad- 
dition to systematic astrometric uncertainties, we es- 
timate a statistical positional error for HD 32297 of 
AO = (4/ 7 r) 1 / 4 /V8hi2 ■ 6 heaxa /SNR ~ 0.06", in agree- 
ment with the formal errors derived from fitting the vis- 
ibilities. Combining the systematic and statistical er- 
rors, the total rms positional uncertainty for HD 32297 
is 0.09". Thus, the peak continuum emission is offset 
from the stellar position at the 4cr level. 



Millimeter Imaging of HD 32297 
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Fig. 3.— Two-population SED fit (see text and F07). The 
small-grain population is plotted in red, the large-grain population 
in blue, the stellar photosphere in grey, and the composite model 
in green. The relative shadings correspond to the 2-d histograms 
of the component SEDs for each model state in the Markov chains. 
Photometry references are as follows: MICHE LLE (F0 7), T-ReCS 
(extr acted from data in lMoerchen et ail l|2007T )'). IRAS l|Mo6r et ail 
[2003) . MIPS (this paper), CARMA (this paper). The model pro- 
vides a good fit to the mid- and far-infrared data but underesti- 
mates the observed mm flux, potentially providing evidence for a 
third grain population. 

As previously noted, the Figure [1] map suggests the 
source is marginally resolved. To quantify this effect, 
in Figure [2l we show the visibility amplitudes for the 
combined data set as a function of projected baseline 
length from phase center. The data arc binned by base- 
line length in three circularly symmetric annular bins, 
with widths chosen to provide approximately the same 
number of visibilities per bin. The errors in amplitude 
represent the standard deviation in the mean of the visi- 
bilities in each bin. Figure [5] shows that the source ampli- 
tudes decline with baseline length, with a deviation from 
a point-source model of approximately 2a. Thus, the 
source appears to be indeed resolved, though additional 
data are needed to confirm this result. In particular, we 
note that uncorrected atmospheric phase errors at mm 
wavelengths can artificially enlarge targets. Still, the ob- 
served morphology agrees with that observed at other 
wavelengths. Fitting an elliptical Gaussian to the visi- 
bilities yields a FWHM of 1.8±0.6" x 0.2±0.6" with a po- 
sition angle of 55 ±10°, consistent to within la of the ori- 
entation observed in the NIR and MIR (P.A.= 45 - 57°, 
Schneider et al. 2005, F07, Moerchen et al. 2007). 

We observe no CO J = 2— 1 line emission in our 
data. The 3a limit is 22 mJy beam -1 (1.5 K, 6*bcam = 
0.70" x 0.52") i n a 42 km s" 1 chan nel. Adopting the same 
assumptions as lDent et al.l (|2005f ). the corresponding gas 
mass limit is M gas < 0.3Aij up . 

4. DISCUSSION 

F07 attempted to model the observed SED and N'- 
band image of HD 32297 with a single ring of grains of 
characteristic size, but found that a second population of 
grains was needed to adequately fit the observed SED for 
A > 25/im. Indeed, the contrast between the observed 
mm and iV'-band morphologies (see below and Figure [4]) 
suggests that the grains responsible for the emission at 
each wavelength constitute separate populations. To test 
whether the population of mm-emitting grains is consis- 
tent with the second, larger grain population proposed 
by F07 to fit the SED at 25 /mi < A < 60 /im, we revisit 
the F07 model, adopting their data and fitting method, 
and incorporating the Qa-band flux of iMoerchen et aT] 



(2007) and the new MIPS and CARMA fluxes. The free 
parameters in the single population model of F07 are 
the disk inclination (i), position angle (PA), radii of the 
inner and outer edges (wq, w\), surface density power- 
law index (7) , vertical optical depth to absorption at the 
inner edge (r = T^ bs (n7 )), stellar flux factor (£), and ef- 
fective grain size (A sm ). To fit the long- wavelength SED 
(A > 25/im), we augment this model with a population of 
larger grains of effective size Ai g and total emitting area 
A\ g , located in a narrow ring at the small-grain inner-disk 
edge, tuq. These grains contribute flux, F u ,\ g , according 
to: 

Aig/Aif A > Ai g , 
1 otherwise, 



(1) 



Tig(r) =468 



A, 



Aig/1 /im 



1/5 



1 AU 



-2/5 



K, 



Fu^=[-g-)e v>]s B v [T^(w )]. 



(2) 
(3) 



Following F07, we ran three Monte Carlo Markov 
chains with 3 x 10 4 samples, and simulataneously fit both 
the large and small grain populations. The results of 
this procedure are listed in Table [H and the correspond- 
ing range of allowed dust emission is plotted in Figure [3] 
(red: small grains, blue: large grains, grey: photosphere, 
green: composite) . While the two-population model pro- 
vides a satisfactory fit to the mid- and far-infrared data, 
the mm flux is underestimated at the 4c level, suggesting 
that the second grain population proposed by F07 is not 
responsible for the majority of the mm flux. Neverthe- 
less, the new MIPS data lend further evidence to the F07 
suggestion that two populations are needed to fit the ob- 
served SED for A < 160/im. Therefore, the fit suggests at 
least three distinct populations are traced by the current 
observations. However, from the two-population model, 
only the 1.3 mm flux appears to trace the putative third 
population. Since this population is described by both 
a mass/emitting area and a size/temperature, we lack 
sufficient data to fully characterize it via modeling. We 
note, though, that this population likely traces > 95% of 
the total dust mass. Assuming a characteristic stellocen- 
tric distance of ~50 AU (§3), L* = 5.4L© (F07), and an 
effective grain size of 1.3- mm, the implied mm grain tem- 
perature is ^30 K, suggesting a dust mass of M mm ~ M e 
(adopting an opacity of 1.7 cm 2 g -1 ). This estimated 
mass is among the highest observed for debris disks de- 
tected in the (sub)mm and is two orders of magnitude 
larger than that implied for the "large-grain" population 
in the SED fit: M Jg ~ 0.02M e (using the fitted param- 
eters in Table [1] and assuming spherical grains with a 
density of 1 g cm -3 ). Future far-infrared / sub-mm ob- 
servations are needed to confirm the three populations 
proposed here and better constrain their properties. 

To address the observed mm morphology, in Fig- 
ure [U we qualitatively compare the CARMA mm map 
(contours, left panel) to the MIR image of F07 (con- 
tours, right pan e l) and the NIR scattered-light image of 
iSchneider et alj (|2005f ) (color, both panels). The images 
at all three wavelengths are consistent with an edge-on 
disk. However, while both the NIR and CARMA data 
exhibit a brightness asymmetry between the northeast 
and southwest lobes inward of ~0.5", the F07 MIR im- 
age is consistent with azimuthal symmetry. We note that 
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TABLE 1 
Best-Fit Model Parameters 





Parameter 


r 

Best-fit 


; — ; 

Description 


i (deg) 


90 ±5 


disk inclination 


PA (deg) 


46 ±3 


disk position angle 


zu (AU) 


70 +20 
' u -10 


inner edge 


zsi (AU) 


> 1200 


outer edge 


logio T 


-2 4 +0 - 3 


vertical optical depth at wq 


7 


< -1.63 


surf density power-law index 


df/7.22 x lCT 20 


1.03 ± 0.03 


(R*/d) 2 , stellar flux factor 


log 10 (Asm/l/im) 


_i 4+0-6 


small grain effective size 


logio(Aig/lcm 2 ) 


28.9 ±0.2 


large grain emitting area 


logioC-W 1 / 1111 ) 


i q+0.2 


large grain effective size 



Note. — Confidence intervals are 95% for marginal poste- 
rior distributions. Adopted priors are as in F07, except for 7, 
which is constrained to be in the domain [-4,0] due to the incon- 
sistency of rising surface density with the scattered-light image 
and CARMA map; the large grain effective size was constrained 
by a log-uniform prior from 1 nm to 100 mm and a requirement 
that Ai g > A sm . 




Fig. 4.— Left: CARMA contours of HD 32297, overlaid on the 
NIR scattered-light image from Schneider et al. (2005). Right: 
Photosphere-subtracted MIR contours from F07 overlaid on the 
same image. The asymmetry in the CARMA data between the 
northeast and southwest lobes suggests the large, mm-sized grains 
may be trapped in resonance with an unseen exoplanct. A similar 
asymmetry is also observed in scattered light but not in the MIR. 
The contrast between asymmetry at short and long wavelengths 
and symmetry at intermediate wavelengths is a direct prediction 
of the resonant debris disk models of Wyatt (2006). 

iMoerchen et alJ (|2007l ) found evidence for asymmetry in 
their Qa-band data, though their data were not PSF sub- 
tracted, and the asymmetry was in the opposite sense as 
observed in the NIR / CARMA data (NE lobe brighter 

than SW). ^ 

Recently, iGrigorieva et al.l ((2007) used predictions 
from their numerical model of collisional avalanches to 
suggest that the observed scattered-light asymmetry in 
HD 32297 results from the breakup of a large planetesi- 
mal. However, while a massive collision can explain the 
observed NIR and mm morphology, it is not clear why 
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the MIR image would not show a similar asymmetry An 
alternative hypothesis is that the structure results from 
a planetary-induced resonance. In this case, the dust is 
either trapped in resonance as it drifts inward due to 
Poynting-Robertson drag, or it remains locked in reso- 
nance after being generated by parent planetesimals in 
resonance themselves (e.g., Krivov et al. 2007, and refer- 
ences therein). Interestingly , a recent study of the latter 
mechanism by IWvattJ (|2006l ) directly predicts a contrast 
between asymmetry at short and long wavelengths and 
symmetry at intermediate wavelengths, as observed in 
HD 32297. In this model, (sub)mm emission is domi- 
nated by large grains, which have the same clumpy res- 
onant distribution as the parent planetesimals. Small 
grains, traced at short wavelengths, exhibit a similar 
asymmetry, as they are preferentially born in the high 
density, resonant structures before being rapidly expelled 
from the system. Lastly, moderately-sized grains sam- 
pled at intermediate wavelengths remain bound to the 
star, but have fallen out of resonance due to radiation 
pressure and are subsequently scattered into an axisym- 
metric morphology. 

This proposed scenario of IWvattJ (|2006l) conveniently 
explains the qualitative picture for HD 32297 depicted 
in Figure [4l yet we caution that this suggestion is highly 
speculative, and future rigorous modeling of this sys- 
tem is needed to draw firm conclusions from the avail- 
able data. One potential problem with this hypothesis is 
that the SED model predicts that the small (TV'-band- 
emitting) grains have sizes < 1 /im (Table [1]), similar to 
that expected to produce the NIR scattered-light image 
(see discussion in F07). However, the interpretation of 
Figure H in terms of the iWvattl (|2OO60 models requires 
that the NIR-scattering and MIR-emitting grains have 
different sizes. This ambiguity could potentially be re- 
solved through simultaneous modeling of scattering and 
emissi o n, inc orporating the optical image presented in 
IK alas! (|2005[ ). Most importantly, and independent of 
speculation, HD 32297 is currently one of only a few de- 
bris disks with resolved observations in four wavelength 
regimes (optical, NIR, MIR, mm). Taken together, these 
observations can provide a unique testbed for theories of 
grain and planetary dynamics. 
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